arXiv:hep-ex/9701012vl 22 Jan 1997 


CLNS 96/1450 
CLEO 96-21 


Tau neutrino helicity from 


energy correlations 


CLEO Collaboration 

(February 7, 2008) 


Abstract 


We report a measurement of the magnitude of the tau neutrino helicity from 
tau-pair events taken with the CLEO detector at the CESR electron-positron 
storage ring. Events in which each tau undergoes the decay r—with h a 
charged pion or kaon, are analyzed for energy correlations between the daugh¬ 
ter hadrons, yielding \^h\ = 1-03 ± 0.06 ± 0.04, with the first error statistical 
and the second systematic. 
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In the Standard Model of the weak interaction the tan is a seqnential lepton with a pnrely 
left-handed weak isodonblet partner, i/,-; i.e., the tan nentrino has its spin anti-parallel to 
its momentnm. For the decay of the tan to a single psendoscalar hadron h (with h = tt or 
K), the parent tan spin is maximally correlated with the danghter hadron direction |I|. At 
high beam energies, the s-channel interaction > 7 *—> tends to have 

the tan spins aligned, leading to a correlation between the angles of the hadron momentnm 
vectors as measured in the tan rest frames. The Lorentz boost to the laboratory frame 
also maps these angles into the observed particle energies, transforming the spin-direction 
correlation into an energy-energy correlation between the hadrons. In this Brief Report we 
describe a measurement of the magnitude of the tan neutrino helicity based on psendoscalar 
energy-energy correlations at a center-of-mass energy of Y^=2F^b~10.6 GeV. This study 
complements other recent measurements that involve other modes, other techniques, and/or 
higher beam energies p|-§]. 

The helicity, denoted here as and elsewhere as is related to the parameter jav and 
the charged current couplings Qa and gy by 

— hu^ = 'yav = 2gagv/{g^^ -I- gy) (1) 


If we define c± as the cosine of the angle of the momentum in the parent rest frame 
with respect to the boost direction, the correlation is given by P 




a 


dc^dc^ 




( 2 ) 


( 2^2 + ^ 2 ) 

At -^^=10 .6 GeV, the net spin alignment is roughly 90%. Using instead the laboratory 
variables x± = ^Jp% + /Eh gives the energy-energy correlation 


d‘^a 


dx^dx- 


= Fi{x+,X-)+^l ■ F2 {x+,X-), 


( 3 ) 


with Fi 2 being known kinematic functions. This correlation is maximal and identical for 
a left-handed 1 ) and a right-handed (^/i=-|-l) tan neutrino, and it vanishes for no 

preferred handedness (^h=0). It is not explicitly parity-violating but is a consequence of 
parity violation in tan decay. 


The data analyzed here correspond to a luminosity of 1.64 fb ^ (1.5x10® 


r' r 


events), 


collected at the Gornell Electron Storage Ring (GESR) with the GLEO II detector l]!^. A 
set of three concentric drift chambers in a 1.5 T axial magnetic held measures charged par¬ 
ticle momenta with resolution ap/p (%) ~ Y^^((hl5p)^'^|-~((h5)^, p in GeV/c. Surrounding the 
drift chambers, but inside the superconducting magnet coil, is a GsI(Tl) crystal electromag¬ 
netic calorimeter. Barrel crystals surround the tracking chambers, covering | cos 6^1 <0.82, 
with 6 the angle with respect to the e~^ beam direction. Two identical endcaps occupy 
0.80<| cos6'|<0.98. Particle time-of-hight is provided by 5 cm-thick scintillation counters 
located just inside the calorimeter in the barrel and endcap. Muons are identihed by their 
penetration through the calorimeter, coil and one or more of three 36 cm-thick slabs of 
magnet iron; three layers of larocci tube chambers instrument the gap behind each slab. A 


three stage hardware trigger []IT| uses combinations of calorimeter, tracking chamber, and 
time-of-hight information to initiate detector readout. 
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To measure we first extract —>■ {h^T7)[h~u) events that are relatively free of 

backgrounds from non-tan processes. Monte Carlo samples of events from V and V±A 
processes are then generated and mixed appropriately to obtain a desired value of ^h- We 
determine what value of best hts the data using likelihood and methods, investigating 
both two-dimensional (x+,a;_) and one-dimensional (c+c_) distributions. The techniques 
are then analyzed for systematic biases and uncertainties. 

Events are selected that have exactly two oppositely charged tracks in the hducial volume 
of the tracking system with |cos6'|<0.7 and a;±>0.2 and which project back to the 
luminous region. To reduce the contamination of non-rr QED events, we demand that at 
most one track have a;±>0.8, that neither track have a;±>0.95, and that the total shower 
energy in the calorimeter be less than 0.85^^. To remove backgrounds from hnal states that 
involve electrons, each track must have its calorimeter energy less than 85% of its measured 
momentum. These criteria reduce the number of events to about 115,000. 

Final states that contain photons from neutral pions are suppressed by demanding that 
there be no photon-like showers in the calorimeter that are not matched to charged tracks. 
A combination of energy and isolation criteria is used to identify such showers. This require¬ 
ment also greatly reduces the contamination from events with a single photon, such as from 
the process e'^e~^pj. To suppress these radiative events further the two charged tracks are 
required to have an opening angle of greater than 90° in the r — (p projection. To minimize 
possible systematic effects, events must satisfy trigger criteria that are suitable for two-track 
events and that are uniform throughout the dataset. 

The most effective variable for removing events from 77 interactions is Omin, dehned by 


sin© 


min 


\p±\ _ 

Eb ■ {2 — X- — x+) ’ 


(4) 


i.e., the ratio of \p±\, the missing net momentum transverse to the beams, to the missing 
energy of the event 10 . 0 . This requirement also effectively suppresses radiative QED events 
in which radiated photons go undetected. Events are retained if sin0min>O.lO. Because 
particles from 77 interactions tend to have low momentum, at least one of the tracks is also 
required to have a;±>0.3. 

Discrimination against muons involves four mutually exclusive criteria, three of which 
use information from the muon chambers and one of which demands the particle leave an 
energy deposition in the calorimeter that is inconsistent in magnitude and shape with that of 
a muon. This last criterion is invoked only if the charged track projects into uninstrumented 
regions at the azimuthal boundaries of the iron absorber. Each track must pass one of these 
four criteria for the event to be accepted. The efficiency of the combination of these criteria 
to select a track as a pion has been measured using an independent subsample of hadrons 
in T^T~ decays involving a lepton recoiling against one or more neutral pions and a single 
track that is assumed to be a charged pion or kaon [i — hlrni^)). This efficiency rises sharply 
from a threshold at x± ~ 0.2 to 60% at x± ~ 0.25, and reaches a plateau of about 70% by 
x± ~ 0.5. Studies of radiative mu-pair (p/ry) events in both data and simulation indicate a 
rate for misidentihcation of a muon as a pion of less than 2 % per track over the kinematic 
range of this analysis. 

The resulting 2041 h^h~ candidates are binned in x+ vs. X- with bin size 0.1 and range 
0.2<x±<0.95, as shown in Fig. ^ The missing corner bins at x+=x_ are the result of the 
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energy criteria that help suppress 77 events and QED events. This hgure clearly shows the 
depletion in the other corners that have one “hard” and one “soft” hadron, as expected for 

V ± A. In addition we examine the product c+c_ as shown in Fig. ^ 

The Monte Carlo simulations start with the KORALB |T^ tan-pair generator and use the 
GEANT package to model the detector response. Monte Carlo studies indicate that over 
97% of these events were of the five final states 7 r+ 7 r“z/P (81%), (9%), (4%), 

/u) ( 2 %), and 'k^K*^vu ( 2 %); the next leading mode was 7 r=‘=e^z/p(z//p) ( 0 . 8 %). 
For these five specihc decay channels large samples were generated with both iV—A) and pure 

V coupling. In the case of events were also thrown with a iV+A) coupling; 

for the other hnal states the iV+A) and {V—A) predictions are identical. These events are 
then appropriately scaled to be added together for comparison to the data distributions. For 
each bin in either the (a;+, x_) or c+c_ analysis, the number of predicted events for any given 
value of is expressed in terms of the number from the iV—A), {V+A) and V simulations 
by 


+ G+iiVf+ GoiV^ 


( 5 ) 


with C±i=(^l ± ^h)/2 and C'o=(l — A binned maximum likelihood fit in the one pa¬ 
rameter, is then performed, with the resulting likelihood distribution for the (a;+,a;_) 
analysis shown in Fig. i The result is \^h\ = 0.998lo;o59, with the uncertainty being purely 


g, yielding = 


statistical. A similar likelihood curve is found for the c+c_ analysis 

1 ni 1 +0-064 
i.Uii_o,o 60 - 

We have made many systematic checks which can be divided into four broad categories: 
possible biases in the methodology, our estimate of backgrounds, dependence on the simu¬ 
lation of the variables used in event selection, and possible biases in the simulations. These 
systematic effects have been studied for both the (a;+,x_) and c+c_ analyses. Because the 
values of |^/i| from the two analyses are statistically indistinguishable, only the details for 
the {x+,x_) analysis will be given. 

In addition to the likelihood technique, a analysis was performed, with minimal shift 
in the value of Further, studies with very large samples of simulated events (which, for 
speed, had their kinematic variables smeared with typical resolution functions instead of full 
detector simulation) showed that the choice of the likelihood estimator causes no inherent 
shift in the value of Various choices of bin size were tried for both the likelihood and 

analyses, from which we assign a systematic uncertainty of ± 0 . 02 . 

Studies with Monte Carlo samples indicate a bias of —0.011±0.006 in the fitting procedure 
and of —0.020 ± 0.010 due to the use of only the five most common modes. We include these 
shifts to obtain a central value of |^/j| = 1.029 and incorporate the uncertainties on the shifts 
into the overall systematic error. 

If the Monte Carlo were not properly handling the various decays or if there were 

non-trivial amounts of non-tan backgrounds, then changing the tan branching fractions would 
affect the result because the various hnal states have different distributions in the {x+,x-) 
plane. Variations of 10-20% in these branching fractions produced changes in of less than 
0.01. We used an extended likelihood technique to incorporate the normalization into the 
ht, which yielded a tan branching fraction and statistical uncertainty of 13,^ = 0.11 ± 0 . 01 , 

IHM 


consistent with the established value 


This extension of the likelihood to include 


normalization had a negligible effect on the htted value of \^h\, changing it by 0 . 01 . 
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The effect of possible contamination from the /i/i 7 hnal state was investigated by intro¬ 
ducing simulated events into the data sample. After weighting events that passed all 

but the muon rejection criteria by roughly five times the rate at which muons fake hadrons, 
a shift in of —0.02 is observed. The resulting uncertainty on assuming the actual 
measured misidentification rate, is less than 0.005. 

The Monte Carlo simulations for the final state 7r~^7T~i>u were examined to see if any 
of the efficiencies of the various selection criteria were different for V and V—A. Of all the 
selection criteria, the restriction on sin ©min is the one to which is most sensitive, especially 
because it tends to eliminate events along the diagonal in the (a;+,x_) plane. We therefore 
varied the limit on sin0min down to 0.05 and up to 0.20, observing changes in the fitted value 
of <0.02. Replacing the limit on sin0min with requirements on acoplanarity or visible event 
energy (also effective in suppressing 77 backgrounds, but less so) produced similar changes 
in but of the opposite sign. An uncertainty of ±0.01 was assigned to this effect. 

Other event selection variables were also evaluated by varying the requirements, leading 
to very small uncertainties. An exception was the h//i discrimination, for which we tried 
using various subsets of the four criteria as well as changing the limits in them individually, 
leading to an associated uncertainty of ±0.025. Similarly we tightened and loosened (where 
possible) the trigger requirements with effects observed of ±0.015. 

All of the efficiencies depend on the Monte Carlo simulations. The most critical aspects of 
this are the photon veto (for which hadronic splitoffs need to be properly modeled), the h/^ 
discrimination, and the trigger. To investigate these more thoroughly, samples of events were 
obtained in both data and in Monte Carlo simulation for the £-h{nn^) topology previously 
described. From these studies we obtained a:±-dependent corrections to the efficiencies for 
all three aspects of the analysis. Applying these indicated very small shifts, some positive 
and some negative, in the fitted value of Based on this we applied no shift in the central 
value and assigned a systematic uncertainty of ±0.01. 

Taking all the effects in quadrature yields an overall systematic uncertainty of = 
±0.040 for the (a;+,x_) analysis. For the c+c_ analysis, the biases shift the central value 
to = 1.036 with a systematic uncertainty of A^h = ±0.049. Our quoted value for the 
helicity and its uncertainties is the simple average of these two analyses. 

In summary, we have measured the magnitude of the helicity of the tan neutrino from 
the energy-energy correlation in r+r“—>(h+I7)(h“i/) events, obtaining |^/i|=1.03±0.06±0.04. 
This is consistent with the other measures of the Ur helicity M- Considering only V and 
A interactions, the physical bounds are 0<|^/i|<l, and our result corresponds to |^/i|>0.87 
at 95%CL. From Eqn. 5 this implies at most 24% pure V (or pure A) at 95% CL. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun¬ 
dation, the U.S. Department of Energy, the Heisenberg Foundation, the Alexander von 
Humboldt Stiftung, Research Corporation, the Natural Sciences and Engineering Research 
Council of Canada, and the A.P. Sloan Foundation. 
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FIG. 1. The energy-energy distribution for the data, uncorrected for efficiencies, using the 
scaled energy variables. Note that the actual range of the data is 0.2<x±<0.95. 
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FIG. 2. The distributions of c+c_ for data and for V—A and V simulations. 
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FIG. 3. Log likelihood vs. for the (x+,x_) analysis. The maximum log likelihood value 
has been subtracted. The insert shows the full range explored. 
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